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Abstract. The evaluation of quantum corrections in the theory of the electroweak and strong interactions via higher- 
order Feynman diagrams requires complicated and laborious calculations, which however can be structured in a strictly 
algorithmic way. These calculations are ideally suited for the application of computer algebra systems, and computer 
algebra has proven to be a very valuable tool in this field already over several decades. It is sketched how computer 
algebra is presently applied in evaluating the predictions of the electroweak theory with high precision, and some recent 
results obtained in this way are summarized. 



INTRODUCTION 

The electroweak and strong interactions of elementary 
particles are very successfully described by quantized 
gauge field theories. The quantized nature of these theo- 
ries manifests itself via corrections beyond the lowest or- 
der in the perturbative expansion, which is based on Feyn- 
man diagrams. The evaluation of higher-order Feynman 
diagrams (which are called loop diagrams) is a techni- 
cally very complicated but on the other hand algorithmic 
procedure. The development of computer-algebra sys- 
tems was boosted by the demand for this kind of applica- 
tions in particle physics, and Schoonschip (1) was one of 
the first implementations of a powerful computer-algebra 
program. Further examples of computer-algebra systems 
that have their roots in particle physics are Reduce (2), 
Macsyma (3), Mathematica (4), and FORM (5). 

Computer-algebra systems allow to perform symbolic 
manipulations and algebraic calculations without round- 
off errors. They are equipped with a number of built-in 
algorithms which provide the basis for the user to imple- 
ment his/her own algorithms for handling specific prob- 
lems. Present-day computer-algebra systems furthermore 
possess capabilities for communicating with external pro- 
grams, e.g. with routines for numerical evaluation, text 
processing and graphics, or with other computer-algebra 
programs. 

Examples of computer-algebra systems being widely 
used at present in high-energy theory are Mathemat- 



ica (4), Maple (6) and FORM (5). Mathematica and 
Maple are general-purpose programs containing a large 
number of built-in functions (and many additional soft- 
ware packages are available). These programs offer ca- 
pabilities for both symbolic and numerical computations, 
support graphical display, and possess user-friendly in- 
teractive platforms. The application of these systems to 
problems in high-energy physics involving expressions 
with a huge number of terms can be limited by the com- 
puting speed or by memory problems. The latter applies 
in particular to non-local operations, like e.g. factoriza- 
tion, which require to have all terms of an expression 
available within the physical memory of the computer. 
FORM, on the other hand, is a program that was specif- 
ically optimized for handling very large expressions. It 
is less user-friendly than Mathematica and Maple, con- 
taining much fewer built-in operations and allowing only 
non- interactive execution. For recent developments con- 
cerning the parallelization of FORM, see Ref. (7). 



PERTURBATIVE EVALUATION OF 
GAUGE THEORIES WITH 
COMPUTER-ALGEBRAIC METHODS 

The concept of treating interactions as a perturbation 
to a free field theory and performing an expansion in the 
coupling constants leads to a description of scattering 



processes in terms of Feynman diagrams. The lowest- 
order prediction, corresponding to the classical limit, for 
a process with a certain number of external particles (e.g. 
a 2 — > 2 scattering process) is obtained from the sum 
of the connected diagrams containing the lowest pos- 
sible power of the coupling constants (which enter via 
the interaction vertices). These are in general tree-level 
diagrams. In higher-order diagrams additional interac- 
tion vertices give rise to closed loops of propagators, for 
which an integration over the internal momenta has to be 
performed. 

The prediction for a scattering process of certain fields, 
assigned to the external legs, and a specified number of 
loops can be obtained via an algorithmic procedure. In a 
first step, all topologically different diagrams (for which 
in renormalizable theories only 3-point and 4-point in- 
teraction vertices are possible) have to be generated. In- 
serting the fields of the model under consideration into 
the topologies in all possible ways leads to the Feynman 
diagrams. The Feynman rules translate these graphical 
representations into mathematical expressions. 

Since the loop integrals in general lead to divergences, 
the expressions need to be regularized (i.e. made math- 
ematically meaningful). In a renormalizable theory the 
divergences can be absorbed into a redefinition of the pa- 
rameters of the theory. The renormalization is further- 
more necessary in order to fix the physical meaning of 
the parameters order by order. 

The evaluation of the Feynman ampUtudes involves a 
treatment of the Lorentz structure of the amplitude, cal- 
culation of Dirac traces etc. At the one-loop level it is 
possible to reduce all tensor integrals to a set of stan- 
dard scalar integrals, which can be expressed in terms of 
known analytic functions. As a consequence, with the 
existing techniques a wide class of processes with up to 
four external legs can be evaluated at the one-loop order 
in massive gauge theories (for a discussion of the tech- 
nical problems occurring in one-loop processes with six 
external legs, see e.g. Ref. (8)). 

In contrast to the one-loop case, no general algorithm 
exists so far for the evaluation of two-loop corrections in 
the electroweak theory. The main obstacle in two-loop 
calculations in massive gauge theories is the complicated 
structure of the two-loop integrals, which makes both the 
tensor integral reduction and the evaluation of scalar inte- 
grals very difficult. In general the occurring integrals are 
not expressible in terms of polylogarithmic functions. For 
the evaluation of some types of integrals that do not per- 
mit an analytic solution numerical methods and expan- 
sions in their kinematical variables have been developed. 

Applying the appropriate on-shell conditions to the ex- 
ternal legs one obtains the S-matrix element from the sum 
of all contributing Feynman amplitodes. Squaring it and 



performing the phase space integrations one finally ar- 
rives at predictions for cross sections and life times. 

Computer-algebraic methods can faciUtate most of the 
above-mentioned steps. Besides benefits from automa- 
tion, a computer-algebraic treatment is also useful for ver- 
ifying the correctness of the different steps of a certain 
calculation. In particular, results obtained at the algebraic 
level (before inserting specific numerical values for the 
parameters) are well suited for highly non-trivial checks, 
e.g. with respect to their UV- and IR-finiteness, gauge- 
parameter independence, and the validity of Slavnov- 
Taylor identities. As an example, in Ref. (9) a Slavnov- 
Taylor for the two-loop Z-boson self-energy in the elec- 
troweak Standard Model (SM) has been verified by show- 
ing that the results of about 4000 Feynman diagrams add 
up to zero algebraically. 

As indicated by the above example, powerful 
computer-algebraic tools are very useful for calculations 
(in particular of higher-order corrections) in the SM, since 
the large number of different fields in the SM gives rise 
to a large number of contributing Feynman diagrams (at 
the one-loop level typically (10^), at the two-loop level 
0(10^)), and the massiveness of the fields makes the eval- 
uation of the loop diagrams very compUcated in general. 
The technical complications are even higher in exten- 
sions of the SM. In the Minimal Supersymmetric Stan- 
dard Model (MSSM) the duphcation of the number of 
fields compared to the SM leads to a plethora of possible 
interaction vertices and consequently to a large increase 
in the number of diagrams contributing at a certain or- 
der. In QCD, on the other hand, computer-algebraic tools 
are particularly valuable for multi-loop applications. In 
Ref. (10), for instance, the four-loop P function of QCD 
has been calculated. This required the computation of 
about 50000 diagrams, showing clearly the need for a 
high degree of automation. Similarly, for tree-level pro- 
cesses with many particles in the final state thousands of 
diagrams can contribute and algebraic methods can be 
useful for obtaining compact and numerically efficient 
representations, see e.g. Ref. (11). 

Examples of computer-algebra based collections of 
program packages presently used for higher-order calcu- 
lations in the electroweak theory and QCD are (where the 
different programs in each collection mostly use common 
syntax and can be linked together) 

(i) PeynArts (12), PeynCalc (13), FormCalc (14), 
TwoCalc (9, 15), LoopTools (14), s21se (16), 

(ii) GEFICOM (17), QGRAF (18), MATAD (19), 
MINCER (20), 

(iii) DIANA (21), QGRAF (18), ON-SHELL2 (22), 

(iv) xloops (23), GiNaC (24). 



FeynArts, FeynCalc, FormCalc and TwoCalc are writ- 
ten in Mathematica {FormCalc is partially written in 
FORM). FeynArts is a program for generating all Feyn- 
man amplitudes contributing to a certain process to a 
given order in Mathematica format and for drawing the 
corresponding Feyrmian diagrams. As a feature of partic- 
ular importance for higher-order calculations in the elec- 
troweak theory, FeynArts generates not only the unrenor- 
malized diagrams at a given order but also the coun- 
terterm contributions at this order and the counterterm 
diagrams needed for the subloop renormalization. The 
model files for the electroweak SM (including the Feyn- 
man rules for the background-field formulation of the 
SM (25)) and QCD are predefined in FeynArts. Re- 
cently also the model file for the MSSM has been com- 
pleted (26). In applications to other models, e.g. chiral 
perturbation theory (27), the appropriate model file has to 
be provided by the user. This was also the case in previ- 
ous applications in the two Higgs-doublet model (28) and 
the MSSM (29,31,32, 33). 

FeynCalc and FormCalc are programs (using the 
FeynArts syntax) for algebraically evaluating one-loop 
diagrams in the electroweak theory and QCD with up to 
four external legs in a highly automatized way. Form- 
Calc internally uses an interface to FORM, which is used 
for the memory- and time-intensive parts of the calcu- 
lation. FormCalc can directly be linked to LoopTools, 
which contains routines for the numerical evaluation of 
scalar one-loop integrals and one-loop tensor coefficients. 
LoopTools is based on the FF (34) package and provides 
a Fortran and a C++ hbrary. 

As mentioned above, much less tools are available 
for two-loop calculations in massive gauge theories com- 
pared to the one-loop case. The program TwoCalc is 
based on an algorithm for the tensor reduction of gen- 
eral two-loop 2-point functions, which extends the algo- 
rithm for the tensor reduction of one-loop integrals (35). 
TwoCalc can be used for an automatic reduction of Feyn- 
man ampUtudes for two-loop self-energies with arbitrary 
masses, external momenta, and gauge parameters to a 
set of standard scalar integrals. It can directly be linked 
to the program s21se, which is written in C++ and per- 
forms the evaluation of the scalar two-loop 2-point inte- 
grals by means of one-dimensional integral representa- 
tions in terms of elementary functions, which allow a fast 
numerical evaluation with high precision. 

The program collection GEFICOM, QGRAF, 
MATAD, and MINCER is mainly used for calculations 
in QCD and for the evaluation of QCD corrections to 
electroweak observables. GEFICOM acts as the master 
program that calls the other packages. It contains Mathe- 
matica and Fortran routines as well as elements written 
in the script languages AWK and PERL. For details of 



GEFICOM and some examples of its applications, see 
Ref. (17). 

The Fortran program QGRAF is an efficient genera- 
tor for Feynman diagrams. As output the diagrams are 
encoded in a symbolic notation. Being optimized for 
high speed, QGRAF is particularly useful for applica- 
tions involving a very large number (i.e. o(10'*)) of di- 
agrams. Within GEFICOM, the evaluation of the dia- 
grams proceeds by performing expansions in their kine- 
matical variables. The resulting integrals are then com- 
puted with MINCER and MATAD. The program MIN- 
CER performs the computation of integrals up to three- 
loop order where all lines are massless and only one ex- 
ternal momentum is non-zero. It makes in particular use 
of integration-by-parts methods (36). While its original 
version was written in Schoonschip, the present version 
of MINCER is reaUsed in FORM. The FORM program 
MATAD was designed for the computation of vacuum 
integrals up to three-loop order which contain only one 
mass scale (i.e. their propagators are either massless or 
carry a common mass). 

The C program DIANA is designed as a master pro- 
gram for higher-order calculations, i.e. it calls the nec- 
essary subprograms for a specific computation. It reads 
the output of QGRAF and can produce a graphical rep- 
resentation for the diagrams if the relevant topologies are 
pre-defined by the user. For the calculation of the dia- 
grams FORM programs are called, e.g. the package ON- 
SHELL2 which can be used for the calculation of single- 
scale two-loop 2-point functions (diagrams with only one 
non-zero mass in the internal lines and the external mo- 
mentum on the same mass shell). 

xloops is a Maple package for calculating certain 
one-loop and two-loop diagrams in the electroweak SM, 
which is linked to C++ routines for numerical integration 
of loop integrals. The symbolic part of xloops is planned 
to be based in the future on GiNaC, which is a specifically 
designed framework written in C++. 



EXAMPLES OF HIGHER-ORDER 
RESULTS IN THE SM AND THE MSSM 

In the following some examples are sketched of recent 
higher-order results obtained with FeynArts, TwoCalc 
and s21se in the electroweak SM and the MSSM. Within 
the SM, higher-order calculations are necessary for the 
comparison of the theory predictions with the experimen- 
tal results for electroweak precision observables like Mw, 
sin^9eff etc. which have meanwhile reached an accuracy 
of better than 1 x 10~^ (37). The precision tests of the 
SM allow in particular to set constraints on the mass of 
the Higgs boson, which is the last missing ingredient of 



the SM and plays a crucial role for a consistent descrip- 
tion of massive particles. 

In Ref. (39) the currently most accurate prediction for 
the W-boson mass, Mw, within the SM has been obtained. 
It contains it particular the complete fermionic contribu- 
tions at the two-loop level, which are treated exactly, i.e. 
without an expansion in the top-quark or the Higgs-boson 
mass. The result for Mw is shown in Fig. 1 as a func- 
tion of the Higgs-boson mass, Mr. It is compared with 
the current experimental value for Mw (37). The present 
95% C.L. lower bound on Mh from the direct search at 
LEP of Mh =114.1 GeV (38) is also indicated. The plot 
shows the well-known preference for a light Higgs bo- 
son within the SM. Confronting the theoretical prediction 
(allowing a variation of the top-quark mass, nii, which at 
present dominates the theoretical uncertainty, within la) 
with the la region of M^P and the 95% C.L. lower bound 
on Mh, one finds that at the present level of accuracies the 
la regions do no longer overlap. 
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FIGURE 1. The prediction for Mw as a function of Mh for 
mi = 174.3 ±5.1 GeV is compared with the current experi- 
mental value, M^P = 80.451 ±0.033 GeV (37), and the ex- 
perimental 95% C.L. lower bound on the Higgs-boson mass, 
Mh = 114.1 GeV (38). 



By comparing the SM predictions for the precision ob- 
servables with those of extended models, it can be investi- 
gated whether the data allow a distinction between differ- 
ent kinds of possible models. In Fig. 2 the predictions for 
Mw in the SM and the MSSM are shown as a function of 
nit. The MSSM prediction contains the dominant SUSY 
contributions of o(aas) and 0{a?) to the p parame- 
ter (29, 30). The allowed region in the SM corresponds 
to varying Mh in the interval 1 14 GeV < Mh < 400 GeV, 
while in the region of the MSSM prediction the SUSY 
parameters are varied, taking into account the constraints 
from direct searches for SUSY particles. As indicated 
in the figure, the predictions in the SM and the MSSM 



give rise to two bands with only a relatively small overlap 
region. This region corresponds to the SM with a light 
Higgs, Mh < 130 GeV, and to the MSSM with heavy 
superpartners, whose virtual contributions decouple from 
the electroweak precision observables. 

The predictions forMw in the SM and the MSSM are 
confronted in Fig. 2 with the current experimental accu- 
racies of Mw and mt (LEP2/Tevatron, outermost ellipse) 
and with the prospective accuracies at the LHC and a fu- 
ture Linear Collider (LHC/LC) and at a high-luminosity 
Linear Collider running in a low-energy mode on the Z- 
boson resonance and the W-pair threshold (GigaZ). As 
can be read off from the figure, the data on Mw and mt 
presently show a slight preference for the MSSM over the 
SM, which however statistically is not very significant. 
The figure shows that the next generation of colliders, in 
particular a Linear Collider in the GigaZ mode, promises 
an enormous improvement in the experimental accuracies 
of Mw and nit (and furthermore also for sin^ 9eff) which 
will allow to test the electroweak theory with unprece- 
dented sensitivity (40). 
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FIGURE 2. The theoretical prediction for Mw within the SM 
and the MSSM in comparison with the current experimental ac- 
curacies (LEP2/Tevatron) and the prospective accuracies at the 
LHC/LC and at GigaZ. 

Besides the indirect constraints from electroweak pre- 
cision tests, supersymmetric models provide a very strin- 
gent direct test since they predict the existence of a rel- 
atively light Higgs boson, whose mass can be calculated 
from the other parameters of the model. 

In Ref. (31) a Feynman-diagrammatic result has been 
obtained for the dominant two-loop contributions to the 
masses of the neutral C 2" -even Higgs bosons in the 
MSSM. The algebraic result obtained with FeynArts and 
TwoCalc has been converted into Fortran code and has 
been implemented into the program FeynHiggs (41). 



While at the tree-level the Ughtest CP -even Higgs 
boson in the MSSM is bounded to be lighter than 
the Z-boson mass, this bound is shifted upwards to 
mh ^ 135 GeV taking into account corrections up to the 
two-loop order. The highest possible values for mh are 
obtained for large values of tanp, the ratio of the vac- 
uum expectation values of the two Higgs doublets of the 
MSSM, large values of the mass of the CfP-odd Higgs bo- 
son, and a large mixing between the superpartners of the 
top quark. Comparing the theoretical prediction for the 
upper bound on mh as a function of tan P with the exper- 
imental exclusion limits obtained at LEP2, it is possible 
to derive constraints on tan p. This is shown in Fig. 3, 
where the excluded region results from combining the 
data of the four LEP experiments (42), and the upper (and 
lower) bound within the MSSM (indicating the boundary 
to the "theoretically inaccessible" region) has been ob- 
tained with FeynHiggs. 

The upper plot shows the case of the so-called "m™'" 
benchmark scenario" (43), in which the MSSM parame- 
ters (for fixed values of mt = 174.3 GeV and the SUSY 
scale MsusY — 1 TeV) are chosen such that mh as a 
function of tanP takes its maximal values. From the 
intersection of the experimentally excluded region with 
the boundary to the theoretically inaccessible region one 
finds an excluded region of 0.5 < tanP < 2.4 within 
this scenario. In the lower plot the MSSM parameters 
have been chosen according to the "no-mixing bench- 
mark scenario" (43), which differs from the m™^ sce- 
nario in that vanishing mixing in the scalar top sector 
has been assumed. In this case a much wider region of 
tanP values, up to about tan P « 10, can be excluded for 
OTt = 174.3 GeV and Msusy = 1 TeV. 



CONCLUSIONS AND OUTLOOK 
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FIGURE 3. The 95% C.L. bounds on mh in the m[J'^'' and the 
no-mixing benchmark scenarios obtained from combining the 
data of the four LEP experiments are compared with the upper 
bound on mh within the MSSM (42). 



The development of powerful computer-algebra sys- 
tems was triggered by applications in high-energy 
physics. Computer-algebra tools have extensively been 
used in this field already for several decades, and many 
of todays calculations would not have been feasible with- 
out computer algebra. A brief overview of computer- 
algebraic methods for the perturbative evaluation of 
gauge theories has been given, and some examples have 
been discussed of recent higher-order results obtained in 
the electroweak Standard Model and its Minimal Super- 
symmetric extension. 

The use of modern computer-algebra programs goes 
beyond their application as tools for certain steps of the 
calculations. As indicated by the above examples of dif- 
ferent collections of programs for precision calculations 
within the theory of the electroweak and strong inter- 



actions, an efficient communication of computer-algebra 
systems with other program components is of particular 
importance. These external programs can be packages for 
numerical evaluations, text processing tools, data bases, 
expert systems, but also other computer-algebra programs 
being particularly well suited for certain sub-parts of the 
problem. In order to facilitate this kind of communica- 
tion, the need for a certain degree of standardization for 
the integration of program parts and the data transfer be- 
tween different systems will become more pronounced in 
the future. 

Accordingly, future improvements of soft- and hard- 
ware promise a further extension of the applicabiUty of 
computer-algebra systems in two different ways. On the 
one hand they could allow highly sophisticated calcula- 



tions which go beyond the scope of present capabiUties. 
On the other hand computer-algebra systems could more 
and more become parts of general problem- solving envi- 
ronments, where different components are integrated in 
such a way that the different parts of calculations and the 
tasks of text processing, graphical representation etc. are 
handled in the most efficient way. 
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